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Abstract 

Various mechanisms have been put forward for cuprate superconductivity, which fit largely into two camps: spin- 
fluctuation and electron-phonon (el-ph) mechanisms. However, in spite of a large effort, electron-phonon interactions 
are not fully understood away from clearly defined limits. To this end, we use a numerically exact algorithm to simulate 
the binding of bipolarons. We present the results of a continuous-time quantum Monte-Carlo (CTQMC) algorithm 
on a tight-binding lattice, for bipolarons with arbitrary interaction range in the presence of strong coulomb repulsion. 
The algorithm is sufficiently efficient that we can discuss properties of bipolarons with various pairing symmetries. 
We investigate the effective mass and binding energies of singlet and triplet real-space bipolarons for the first time, 
and discuss the extensions necessary to investigate d-symmetric pairs. 
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1. Introduction 

The discovery of superconductivity at anoma- 
lously high temperatures in cuprate materials 
caused a major revitalisation in the twin fields of 
correlated electrons and superconductivity. Such 
large transition temperatures could not be ex- 
plained consistently within the BCS approach, 
which is only accurate in the extreme adiabatic, 
weak coupling limit (see e.g. refs [1]). Several mech- 
anisms have been put forward to understand these 
compounds, including very strong el-ph coupling, 
with bipolarons forming a Bose-Einstein conden- 
sate (BEC) in the finite density limit, leading to a 
bipolaronic superconductor (BPS) [2]. 

BPS theory and other approaches such as ex- 
tended Eliashberg theory (EET) [1] are applicable in 
certain limits, either strong coupling (BPS) or weak 
couping (BCS and EET). However, the cuprates 
(and other strongly correlated electron systems) of- 
ten have intermediate coupling between electrons. 
Thus, one cannot expand about a well defined solu- 



tion, and it is necessary to employ advanced simu- 
lation techniques [3,4,5]. Such techniques have the 
power to compute the numerically exact solution 
to a problem, and can provide definitive answers 
about the limits of validity of approximate schemes. 

Here, we simulate a Coulomb-Frohlich (CF) 
Hamiltonian [7] with two electrons of opposite spin, 

H = ~t 4, CT C n(T + ^ V(ll, Vl'^Cnacl^Cn'a 

(nn')<T nn'er 

m m mil " 

Each ion has a displacement £ m . Sites labels are n 
or m for electrons and ions respectively, c annihi- 
late electrons. The phonons are Einstein oscillators 
with frequency ui and mass M. (nn') denote pairs of 
nearest neighbours, and P m = —ihd/d^m ion mo- 
mentum operators. The instantaneous interaction 
V^n, n') has on-site repulsion U and a n.n. attrac- 
tion V (UV model). The Frohlich interaction force 
is long-range, / m (n) = k [(m — n) 2 + l] 3 ^ 2 (k is a 
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Fig. 1. Total energy of the CF model in ID, U/t = 10, 
u/t = 4, V = 0. (Inset: UV model in ID with V = 4. QMC 
results agree with analytic results [8] as a test of the method). 
El-ph coupling is defined as A = ^ m /m(0)/2Maj 2 2:t. 

constant). 

2. Results 

Using the CTQMC algorithm [3], we compute 
the splitting between s and p bipolaron states with 
the estimator E sp = — ln((s})//3. In the p-state, 
s = — 1 for exchanged, 1 for non exchanged con- 
figurations, and s — for paths with cohabiting 
endpoints. Inverse effective mass is calculated as 
mo/m** = {s(Ar) 2 )/(/3h 2 (s)), and Ar is the twist 
in boundary conditions on the imaginary time axis. 
The bipolaron action [5] extends the polaron ac- 
tion [3] . We present triplet splitting in the Coulomb- 
Frohlich model. Below a certain el-ph coupling, the 
singlet-triplet pairing is degenerate, as Coulomb re- 
pulsion stops pairing. The triplet bipolaron is always 
heavier than the singlet state within numerical ac- 
curacy. 

3. Outlook for <i-wave pairing 

We have used an advanced Quantum-Monte- 
Carlo algorithm to simulate bipolarons in the 
presence of very strong Coulomb repulsion with 
long-range el-ph interaction. We have used this al- 
gorithm to investigate the total energy of possible 
bipolaron pairing symmetries for the first time. 

Given experimental evidence for a d-component 
in the cuprates, we aim to simulate the d-bipolaron 
as a prelude to multipolaron calculations. Such as 
state has recently been seen in the two dimensional 
Holstein model as a consequence of large phonon 
frequency in an extended Eliashberg theory [1]. For 
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Fig. 2. Bipolaron inverse effective mass of the CF model in 
ID. U/t = 10, V = and u/t = 4. N.B. Triplets are heavier 
than singlets, as in the solution of the UV model. 

the (i-state, the estimator is similar to the p-state, 
but now parallel end configurations contribute sign 
1 and those rotated by 90° sign -1 [6]. New update 
rules are required in the d-state. The bipolaron ac- 
tion for rotated boundary conditions in time is com- 
plicated, and will be discussed later. While it is clear 
from the estimator for E s( j that a simple bipolaron 
pair cannot have pure d-symmetry, when Coulomb 
repulsion is large, d and s-wave paired states can be 
close in energy. Also, multipolaron effects may be 
capable of inducing a d-wave state [9,1]. 
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